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Abstract 

We add a gravitational background lattice to the simplest holographic model of matter at 
finite density and calculate the optical conductivity. With the lattice, the zero frequency delta 
function found in previous calculations (resulting from translation invariance) is broadened 
and the DC conductivity is finite. The optical conductivity exhibits a Drude peak with a 
cross-over to power-law behavior at higher frequencies. Surprisingly, these results bear a 
strong resemblance to the properties of some of the cuprates. 
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1 Introduction 

Over the past few years, it has been shown that various properties of condensed matter systems 
can be reproduced using general relativity [H El El HI E]- This surprising result is motivated 
by the remarkable gauge/gravity duality which relates a theory of gravity to a strongly coupled 
nongravitational theory. The duality is called "holographic" since the nongravitational theory lives 
in a lower dimensional space. 

The holographic approach has been applied to many interesting phenomena, including super- 
conductivity and superfluidity, Fermi surfaces and non-Fermi liquids. However, with a few notable 
exceptions (see in particular [7j), much of this work has omitted one key ingredient of condensed 
matter systems: the lattice. 

The absence of an underlying lattice becomes particularly important in discussions of optical 
conductivity in systems with a finite density of charge carriers. The underlying translational 
invariance means that the charge carriers have nowhere to dissipate their momentum, resulting 
in a zero frequency delta function in the optical conductivity which obscures interesting questions 
such as the temperature dependence of the DC resistivity. Until now, the most common method 
to avoid this was to treat the charge carriers as probes, their momentum absorbed by the large 
number of neutral fields represented by the bulk geometry [HI [H [10] . 

In Section 2 of this paper, we begin to remedy the situation by constructing a simple gravita- 
tional dual of a system with a lattice. We start with the minimal ingredients needed to describe a 

2 + 1 dimensional system at finite density: four dimensional gravity coupled to a Maxwell field. The 
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equilibrium configuration at nonzero temperature is just a charged black hole. We then add the 
lattice by introducing a neutral scalar field with boundary conditions corresponding to a periodic 
source. This spatially varying field backreacts on the metric and gauge field, imprinting the lattice 
on the bulk geometry. We numerically solve the coupled Einstein-Maxwell-scalar field equations 
to find this gravitational crystaQ 

By perturbing the lattice background, we can compute the optical and DC conductivities 
without the need to work in a probe approximation. Normally, a perfect lattice with no impurities 
has infinite conductivity due to Bloch waves. However, our system naturally includes dissipation 
due to the black hole horizon and we see the expected broadening of the zero frequency delta 
function. 

Our results are described in Section 3. They contain several surprises. At low frequency we 
find that both the real and imaginary parts of the optical conductivity are well described by the 
simple Drude form 



where B and C are constants. Rather strikingly, a power-law fall off with this exponent is seen 
in the normal phase of some of the cuprates exhibiting high temperature superconductivity [TT] . 
(The offset C that we find is apparently not present in these materials.) This result is robust 
against changes in the temperature, lattice spacing, and "strength" of our lattice. We do not 
have a deep understanding of why our numerical experiments on a simple gravitational model 
reproduces results seen in real materials. 

At high frequency, the optical conductivity rises to a positive constant. This is different from 
the cuprates, but agrees with previous calculations without the lattice and is simply a property of 
conformal field theories in 2 + 1 dimensions. 

Finally, the presence of the lattice also renders the optical conductivity finite as the frequency 
co — > 0. This allows us to determine the temperature dependence of the DC resistivity, p = l/cr(0). 
We find that the DC resistivity is strongly dependent on the lattice spacing. Very similar behavior 
was seen recently by Hartnoll and Hofman who performed an analysis of an ionic lattice, created 
by a periodic variation in the chemical potential [T2]. They observed a scaling behavior of the 
form p ~ T 2A where A is a complicated function of the lattice spacing. This scaling can be 
traced to the locally critical AdS2 near horizon region of the Reissner-Nordstrom black hole, with 
the exponent A equal to the dimension of the operator dual to the charge density evaluated at 
the lattice wavenumber. Even though we introduce our lattice differently, we find that the DC 
resistivity precisely reproduces the scaling of [12] at low temperatures. 

x To make the problem computationally manageable, we only add the lattice in one dimension and only compute 
conductivity in the direction of the lattice. 





where the constant K and relaxation time r are determined from the data. However, at interme- 
diate frequency the optical conductivity exhibits a cross-over to power-law behavior 





3 



2 A Holographic Lattice 



The minimal ingredients necessary to compute conductivity in a holographic framework are pro- 
vided by Einstein-Maxwell theory in AdS4. To this we add a neutral scalar field $ which we will 
use to source the lattice. We work with the Lagrangian, 



S= 16^ ldX ^ 9 



R +T2- \ F abF ab - 2V a $V a $ - 4V($) 



(2.1) 



where L is the AdS length scale and F = dA. Our choice of potential corresponds to a massive 
scalar field with mass m 2 = —2/L 2 , 

Vm = - T2 . (2.2) 
The equations of motion derived from the action (2.1) take the following form 

G ab = R ab + ^g ab - 2[V $V 6 $ - V(<S>)g ab ] - (F ac F b c - 9 -fF cd F cd } = 0, (2.3a) 

V F° 6 = 0, (2.3b) 
□$ - V'{<&) = 0. (2.3c) 

Throughout the paper we shall only consider solutions that live in the Poincare patch of 
AdS. We parametrize the holographic radial direction by the coordinate z and impose boundary 
conditions which fix a conformal boundary metric at z = to be of the form, 

ds 2 g = -dt 2 + dx 2 + dy 2 . (2.4) 

We will introduce the gravitational lattice background by providing a spatially inhomogeneous 
source for the neutral scalar field. Near the boundary, $ takes the form 

§->z<i> 1 + z 2 <f) 2 + 0(z 3 ) (2.5) 

According to the AdS/CFT correspondence, 0i should be regarded as the source for the dimension 
two operator dual to 0, say O^, while 02 represents the expectation value (0$). A general, 
inhomogeneous static solution is sourced by <f>i(x, y). However, here we will consider solutions that 
preserve translational invariance in the y direction, with the lattice varying only in the x direction. 
We choose the source 4>i to be 

4>\{x) = A cos(k x) . (2.6) 

We will refer to k as the lattice wavenumber and A as its amplitude, ko is related to the lattice 
size I in the usual way, k$ = 2-k/I. In the rest of this section, we describe these gravitational 
lattices in more detail. 



2.1 The Lattice 



Since our background is both static and translationally invariant in the y direction, the solution 
is co-homogeneity two: it depends only on the coordinates x and z. The most general static, 
electrically charged, black hole solution compatible with our symmetries can be written as 

L2 r 

^2 



ds 2 



1 - z)P{z)Q tt dt 2 + 



Qzzdz' 



P{z){l-z) 



+ Qxx(dx + z 2 Q xz dz) 2 + Qyydy 



(2.7a) 
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with 



and 



$ = z 4>(x, z) 



A= (1- z)if)(x,z) dt 



(2.7b) 



(2.7c) 



where Qij, for ij£ {t,x,y,z}, and if) are arbitrary functions of x and z, to be determined by 



solving Eqs. (2.3). The line element shown above is invariant under reparametrizations of x and 
z, i.e. we have not yet fully specified our coordinate system. We shall address this issue later. 



The factor (1 — z) that appears both in the metric and gauge field ensures that Eqs. (2.7) have a 
smooth non-extremal horizon located at z — 1, provided that Qtt{x, 1) = Q zz (x, 1) and Q 
Q xz , ip and <f> are smooth functions at z — 1. Finally, the factor P(z) is chosen to be 



Q 



P(z) 



1 + z + z i_^K 



(2. 



and controls the black hole temperature 



AttL 



A*? 



8ttL 



(2.9) 



Note that if Q tt = Q zz = Q xx = Q yy = 1, Q a 
Reissner-Nordstrom black hole. 



0, if) — fx — [Mi, we recover the familiar planar 



Finding solutions to the system of equations (2.3), together with the ansatz (2.7), is not a well 



defined problem. The reason being that the Einstein equations do not have a definite character, 
in the PDE sense, if a gauge is not chosen. To overcome this difficulty we will use the DeTurck 
method, which was first outlined in [15] . In this method we change the equations of motion by 
introducing suitable "kinetic terms" for Q zz , Q xx and Q xz , and ensure that any solution to this 
new set of equations is a solution to our initial problem, in a specific gauge. 



2.1.1 DeTurck method 

The DeTurck method is based on the so called Einstein-DeTurck equation, which can be obtained 



from Eq. (2.3a), by adding the following new term 



Gab = ^ab — V( a ^b) — 0, 



(2.10) 



where £ a = g [Tf: d (g) — Tf: d (g)] and T(g) is the Levi-Civita connection associated with a reference 
metric g. The reference metric is chosen to be such that it has the same asymptotics and horizon 
structure as g. For the case at hand, we choose g to be given by the line element (2.7a) with 
Qtt = Qzz = Qxx = Qyy = 1 and Q xz = 0. The DeTurck equation can be shown to be elliptic for a 
line element of the form (2.7a) [15] . 

It is easy to show that any solution to G a b = with £ = is a solution to G^ b = 0. However, 
the converse is not necessarily true. In certain circumstances one can show that solutions with 
£ 7^ 0, coined Ricci solitons, cannot exist [14J. Since the above equations are elliptic, they can be 
solved as a boundary value problem for well-posed boundary conditions and the solutions should 
be locally unique. This means that an Einstein solution cannot be arbitrarily close to a soliton 
solution and one should easily be able to distinguish the Einstein solutions of interest from solitons 
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by monitoring £. In our numerical method we will solve Eq. (2.10 ) directly, together with Eq. (2.3b ) 



and Eq. (2.3c), and a posteriori check that £ = 0, to machine precision. This will correspond to a 



solution of Eqs. (2.3) in "generalized" harmonic coordinates defined by £ = 0. 

In order to show that a solution can be found, we need to guarantee that suitable boundary 
conditions are given, and that these are compatible with £ = 0. 

2.1.2 Boundary conditions 

We begin by discussing the boundary conditions close to the conformal boundary. For a solution 



of Eq. (2.10) with £ = 0, one may compute the asymptotic behavior of the field content close to 



the conformal boundary 



Qtt(x,z) 


= l-j<h (x) 2 + 0(z 3 ), 


(2.11a) 


Qzz(x,z) 


Az 2 

= 1 + —faWfaix) + 0(z 3 ) , 


(2.11b) 


Qxz{%i Z) 




(2.11c) 


Qxx{%i 




(2.11d) 


QyyyXi 




(2.11e) 


<f>(x,z) = 


(j)i{x) + z(j) 2 {x) + 0(z 2 ), 


(2.11f) 


ip(x,z) = 


H+[fi-p(x)]z + 0(z 2 ), 


(2.Hg) 



where \i is a constant that denotes the chemical potential, p(x) is the charge density and <f>i(x) 



are the two possible asymptotic decays for the neutral scalar field <p. As explained in (2.6), we 



introduce the lattice by choosing a source <pi(x) with a nontrivial x dependence of the form 

0i (ar) = A cos(k x) . (2.12) 

Notice that although the scalar field has a lattice wavenumber ho, the stress-tensor is quadratic in 
$, ensuring that the lattice will be imprinted on the metric and gauge field components Qij and 
ip with effective wavenumber of 2 ko- 



Using the expansion (2.11) one can read off the boundary conditions at conformal infinity 



which, as expected, turn out to be of the Dirichlet type: 

Q a (x,0) = Q zz (x,0) = Q xx {x,0) = Q yy (x,0) = 1, 
Q xz (x,0) = 0, (f)(x,0) = <pi(x) ip(x,0) = ji. 

At the horizon (z = 1), regularity demands the following expansion 

Qv(x,z) = Qf{x) + (1 - z)Q ( i f(x) +0((1- zf) 
<j,( Xl y) = {o) (x) + (1 - z)^\x) + O ((1 - zf) , 
i{,(x, y) = ^ {0) {x) + (1 - z)^\x) +0((1- z) 2 ) , 



(2-13) 



(2.14a) 
(2.14b) 
(2.14c) 
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Figure 1: Ajy as a function of the number of grid points N. The vertical scale is logarithmic, and 
the data is well fit by an exponential decay: log(Ajy) = —17.5 — 0.22 N. 

where the W terms can be determined in terms of the terms and their tangential derivatives, by 
solving the equations of motion order by order in (1 — z). The boundary conditions are Dirichlet in 
Q u , because the above expansion demands Qu(x, 1) = Q zz (x, 1), and Robin boundary conditions 
for the remaining variables. 

2.1.3 Numerical method 

We have used a standard pseudospectral collocation approximation in z, x and solved the resulting 
non-linear algebraic equations using a Newton-Raphson method. We represent the dependence in 
z of all functions as a series in Chebyshev polynomials and the ^-dependence as a Fourier series, 
so the x-direction is periodically identified. Our integration domain lives on a rectangular grid, 
(x,z)e(0,27r/fco)x(0,l). 

In order to monitor the accuracy and convergence of our method, we have calculated the charge, 
Q, contained in our integration domain, for several resolutions. Q can be obtained by integrating 
p(x) in the relevant range 

p2ir/k 

Q= dxp(x). (2.15) 

Jo 

We denote the number of grid points in x and z by N and compute A N = |1 — Qn/Qn+i\ for 
several values of N. The results are plotted in Fig. [TJ We find exponential convergence with 
increasing number of grid points, as expected for pseudospectral collocation methods. 

Furthermore, in order to ensure that we are converging to an Einstein solution rather than a 
Ricci soliton we monitor £ a £ a . All plots shown in this manuscript have £ a £ a < 1CT 10 . 



log(A w ) = -17.5-0.22 N 



\ 
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Figure 2: On the left we show Q xz and on the right 0, for fc = 2, Aq — 1, p, — 1.4 and T/p = 0.1. 
Note that has an effective wavenumber of 2ko. 

2.1. 4 Lattice background 

As an example, we show in Fig. [2] a typical result of our numerical code for Q xz and 0. Note 
that the deviation in the metric due to the lattice, Q xz , is small even though the amplitude of the 
oscillating source generating the lattice is Aq — 1. It is important that this small metric correction 
is treated exactly and not just to first order, since otherwise the perturbation that we will introduce 
to compute the conductivity would not feel the lattice. 

After solving the equations, we can study the dependence on temperature, lattice size, chemical 
potential and amplitude, by varying pi, k , p and A , respectively. Due to conformal invariance, 
we can work at fixed chemical potential and vary T(pi)/p, ko/p and A /p, meaning that the 
physical moduli space is only three-dimensional. 

Note that although the operator O^, dual to the scalar 0, is relevant, the scalar perturbation 



decreases in the infrared. This is because the lattice source (2.6) is symmetric about = 0, while 
spatial gradients in the scalar field are suppressed as we approach the horizon. Nonetheless, as is 
clear in Fig. |2j the scalar is not constant on the horizon. In contrast, if we worked at T = 0, the 
scalar field does indeed return to = in the far infrared. In this sense, the lattice perturbation 



(2.6) is irrelevant. 

Our lattices have p/k^ ~ ^(1), which can be interpreted as order one charged degrees of 
freedom per lattice site^J Although the lattice is sourced by a neutral scalar field, it nonetheless 
induces a mild spatial variation in the charge density p(x). On the left hand side of Fig. [3] we 
have plotted the charge density for a lattice with A = 1.5, k = 2, p = 1.4 and pi = 2.12. This 
corresponds to a rather cold lattice with T/p ~ 0.043. In order to better see the effects of the 
non- homogeneities caused by our choice of 0i(x), we have decomposed p(x) as a Fourier series, 

2 This is heuristic as we have only added the lattice in one direction. 
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Figure 3: On the left we show the charge density p(x) and on the right the absolute value of the 
non-zero coefficients of its Fourier series. 

with coefficients On the right panel of Fig. [3] we show the resulting Fourier coefficients. Note 
that nonzero Fourier modes occur only for k equal to multiples of 2k = 4. As mentioned earlier, 
this is because the stress tensor is quadratic in the scalar field, so the lattice seen by the metric 
and Maxwell field has twice the lattice wavenumber. For N x grid points along the ^-direction, 
a Fourier grid can only hope to resolve up to |fc max | < N x /2. It is reassuring that for k = 2, 
and with N x = 42, the highest multipoles have a magnitude smaller than 10~ 10 , meaning that 
our numerical code is capturing the relevant physics. Since higher wavenumbers come from higher 
powers of the scalar field in the nonlinear solution, the Fourier series exhibits an exponential decay 
with increasing wavenumber, which can be best seen in the logarithmic scale used on the right 
panel of Fig. [3j 



2.2 Perturbing the Lattice 



The main purpose of this paper is to explore transport properties in the presence of our lattice. 
According to the AdS/CFT dictionary, on the gravity side, this is mapped into the study of 
perturbations about the lattice background. 

Our first task is to write the equations governing generic perturbations of Eqs. (2.3). We denote 
background fields with hats and expand all fields as 



9ab = 9ab + Kb, A a = A a + b a , $ = $ + T] 



(2.16) 



where h a b, b a and r\ should be regarded as small compared to g ab , A a and respectively. Ex- 
panding equations (2.3) to linear order determines how perturbations propagate in the background 
{g, A, $}. The resulting system of PDEs takes the following form 

3 



-Oh. 



ah 



2R acbd h cd + 2R (n c h b)c + 2V (a V c h 



(a U b)c 



<-b)c 



L 2 



h ab +AV {a r,V b) ^+2V'^)r l g ab +2V^)h ab 
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+ 2f,% c - F ac F M h cd - h -f F cd F cd - 9 -^f cd F cd + 9 -^h cd F cp F/ (2.17a) 



Db„ 



4 2 " 2 

R ac b c - V a (V c b c ) - V c h cd F da - h cd V c F da - F cd V c h da = (2.17b) 

Or] - V a h ab V b § - h ab V a V b $ - V"($)r] = 0, (2.17c) 

where / = db, h ab = h ab — hg ab /2 is the so called trace-reversed metric perturbation, and curved 
brackets acting on indices indicate symmetrization. 

The above system of equations is invariant under the following set of linear transformations 



b a ^b a + V, 



h a b — > h ab 



aX 



(2.18) 



and 



(2.19) 



b a ^b a + f3 c VcA a + A d V d (3 a 

h a b h ab + 2V( a A>)> 

where \ is an arbitrary scalar function and a the components of an arbitrary four-dimensional 
vector. 

The first set of transformations reflects the U(l) gauge freedom associated with electromag- 
netism, whereas the second reflects diffeomorphism invariance. We shall fix these gauge freedoms 
by selecting the Lorentz and the de Donder gauge, respectively 



V a h 



ab 



0. 



V%a = 0. 



(2.20) 



A few words about these gauge choices are in order. One can achieve any of these gauges by starting 
with an arbitrary metric and vector perturbations and then choosing \ and /3 conveniently. This 
can always be done, because one can show that these two gauge conditions imply a wave-like 
equation for \ an d (3. 



So far, our discussion about perturbation theory can be applied to any solution of (2.3): we 



now specialize to (2.7). We want to study the changes in the transport properties induced by 



our lattice. Because our background is invariant under a time translation Killing field dt, we can 
Fourier decompose our perturbations, 



h ab (t,x,y,z) = h ab (x,y,z)e 



-iLdt 



b a (t,x,y,z) = b a (x,y,z)e 



-iojt 



V(t,x,y,z) = fj{x,y,z) e 



-iut 



(2.21) 

Given that the lattice extends along the x-direction, and leaves the y-direction unchanged, we will 
also assume translational symmetry along y. This in turn means that we will focus on perturbations 
with vanishing b y , h ty , h zy and h xy , and that the remaining variables do not depend on y. We are 
left with 11 unknown functions, {h tt , h tz , h tx , h zz ,h zx , h xx , h yy ,b t ,b z ,b x ,fj} in two variables, x and 



z. 



It seems that the resulting system of equations governing perturbations about our lattice 
is overdetermined, because we have 15 non-trivial differential equations to solve, namely the 



{tt, tz, tx, zz, zx, xx, yy} components of Eq. (2.17a), the {t,z, x} components of Eq. (2.17b), the 



scalar equation (2.17c), the {t, z, x} components of the first equation in (2.20) and the last equation 



of (2.20). However, due to gauge invariance, the {tt,tz } yy} components of Eq. (2.17a) and the {t} 



component of Eq. (2.17b) are automatically satisfied if the remaining equations are satisfied. We 
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have explicitly checked that this is indeed the case. We are left with a system of 11 PDEs in eleven 



variables which can be solved using the numerical techniques of Sec. 2.1.3 In the next section we 
shall present how to extract the relevant physics from the resulting numerical perturbations, and 
what boundary conditions do we choose in order to the study the lattice transport properties. 

2.2.1 Boundary conditions for the perturbations 

We wish to study the optical and DC conductivity of the dual CFT in the presence of the lattice. 
The DC conductivity can be extracted from the optical conductivity in the cj — > limit so, in this 
section, we focus on how to extract the optical conductivity only. This issue is naturally connected 
with what boundary conditions we impose for h, b and fj. 

The optical conductivity of the boundary CFT is defined by the following 

A couple of comments are in order about this expression. First, it is manifestly invariant under 
the local U(l) of the electromagnetic field, because it is solely written in term of / = db and 
not b itself. Second, it is also gauge invariant under diffeomorphisms, because f zx is zero on the 



background solution (2.7). Note that f zx can be interpreted as a current, whereas f xt is an electric 
field. Third, the conductivity of the boundary field theory will in general be a function of x. Since 
we will impose a homogeneous boundary electric field, we are interested in the homogeneous part 
of the conductivity which we denote as cr(u). This is the quantity we study below. 

In order to proceed, we first need to generate a nonvanishing f zx and f x t on the boundary. This 
can be easily implemented by requiring b x to be a constant on the boundary, which without loss 
of generality we can set to 1. This rescaling is always possible because we are studying a linear 
system of PDEs. A constant b x will generate a time dependent but homogeneous boundary electric 
field. We also do not want to change the boundary chemical potential, so we will set b t = at 
the boundary. Finally, we do not want to change the lattice, so we will demand fj to vanish as we 
approach z = 0. For the remaining variables we demand normalizability at z = 0. These boundary 
conditions induce the following asymptotic decays for b and fj: 

~b t {x,z) = 0{z), b x {x,z) = l+ Jx {x)z + 0{z 2 ), b z (x,z) = 0{z 2 ) and fj(x, z) = 0{z 2 ). 

(2.23) 



If we input these decays into (]2.22|), we get the following simple expression for the optical conduc- 

(2.24) 



tivity in terms of boundary data only 

~, \ 3x(x) 



l UJ 

For completeness, we also present the decays of the relevant metric perturbations as z — > 

h u (x, z) = O(z), h tx {x, z) = 0(z), h tz (x, z) = 0(z), h xx (x, z) = 0(z), 

h xz (x,z) = 0(z 4 ), h zz (x,z) = 0(z), and h yy (x, z) = O(z). (2.25) 

At the horizon, we impose ingoing boundary conditions, since this is required for causal prop- 
agation. The precise conditions on the perturbation are given in the Appendix. 
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Figure 4: The real and imaginary parts of the perturbation of the scalar field are shown for /i = 1.4, 
T/fi — .115, k = 2, A = 1.5. The two figures on the top have oj/T = 0.06, whereas the two 
figures on the bottom have u/T = 0.6. 
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2.2.2 Numerical implementation of the perturbation theory 



Once we know the boundary conditions both at the horizon and asymptotic infinity, we just have 
to change to new variables which are adapted to the numerics. Because we are solving these 
equations using a pseudo-spectral collocation methods on a Chebyshev grid, we better guarantee 
that no non- analytic terms arise. In order to ensure that, we factor out the non-analytic behavior 



of the near horizon expansion (A.l) from the variables that we actually use in the numerics. For 



example, instead of using fj itself, we work with 

q(x,z) = z(l - z 3 )^)rj(x,z), (2.26) 



with -P(l) given in (2.9). Note that at the conformal boundary q(x,z) has a purely Dirichlet 
boundary condition q(x, 0) = and at the horizon, because the field equations fix fj^(z) in terms 
of the coefficients, q(x, z) has a Robin-type boundary condition. It turns out that we can always 
recast the boundary conditions at the conformal boundary and horizon as homogenous Dirichlet 
or homogeneous Robin boundary condition, by multiplying suitable powers of z and (1 — z 6 ) 3p w 
to our original variables. The only exception to this procedure is b x , which has an inhomogenous 
Dirichlet boundary condition at z = 0. Recall that we need b x (x,z) = 1 in order to generate a 
boundary electric field. 

After discretization, our system of PDEs can be written as a linear map of the form 

M ■ x = x b , (2.27) 

where x b includes the nonhomogeneous boundary condition for b x . This equation can then be 
solved by the LinearSolve in-built function in Mathematica. In Fig. [4] we show the typical output 



of our code for q(x,z) defined in Eq. (2.26). 



3 Conductivity 

Before describing our new results, we first review the conductivity in a translationally invariant 
holographic background. For boundary theories with two spatial dimensions, the conductivity is 
dimensionless and at the conformal point, with [i = 0, is known to be a independent of u, reflecting 
an underlying electron- vortex duality [15]. 

In the presence of a chemical potential, /i, the optical conductivity shows more structure []]. 
Both real and imaginary parts are shown by the dashed, black curves in Fig. [5] for a temperature 
T I \x = 0.115. At large frequency, u ^> ji, the conductivity tends towards the constant, real value 
observed at the conformal point. Indeed, this is the expected behaviour for any scale invariant 
theory in two spatial dimensions. At lower frequencies, u < fi, a drop in Re a reveals a depletion 
in the density of charged states. However, for the purposes of our present discussion, the most 
important feature of the conductivity does not show up in numerical plots of Re a: it is a delta- 
function spike at u = 0. The presence of this delta-function can be seen in the plot of Im cr(o>) 
where, via the Kramers-Kronig relation, it reveals itself as a pole, 

K 

Imer(a;) — > — as u — > (3-1) 

UJ 

for some constant K. 
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Figure 5: The optical conductivity, both without the lattice (dashed line) and with the lattice 
(solid line and data points) for \i = 1.4 and temperature T/fi = 0.115. Note that the lattice 
(which has wavenumber ko = 2 and amplitude Aq = 1.5) only changes the low frequency behavior. 
The pole in Im a without the lattice reflects the existence of a u = delta-function in Re a. 

There is nothing mysterious about the presence of this delta-function. It follows solely on the 
grounds of momentum conservation in the boundary theory. If we have a translationally invariant 
state with nonzero charge density, then one can always boost it to obtain a nonzero current with 
zero applied electric field. This results in the infinite DC conductivity. 

The introduction of a background spatial lattice, as described in the previous section, resolves 
this issue. With no translational invariance, there is no momentum conservation and the u = 
delta-function spreads out, revealing its secrets. In this section we describe what was hiding in 
that delta-function. 

The optical conductivity, a(co), in the presence of the lattice is shown by the solid line in Fig. [5j 
At high frequencies, u 3> fi, the optical conductivity in the lattice background remains unchanged 
from the translationally invariant black hole. The interesting physics lies at lower frequencies. 
The dissipative part of the conductivity, Re a, now rises at low u. This is the redistribution of the 
spectral weight of the delta-function. Moreover, the pole in the responsive part of the conductivity, 
Imcr, has now disappeared, with lma(uj) — > 0, as u — > 0, confirming that there is no longer a 
delta-function at zero frequencjj^J We now describe the characteristics of the conductivity in more 
detail. 

3 The resolution of a delta-function into a Drude-like peak has been seen in a somewhat different context in 
conformal fixed points with vanishing charge density |16j . Here the delta- function is resolved by interactions rather 
than breaking of translational symmetry, either in an e expansion [16] or a 1 /N expansion |17j . 
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Figure 6: A blow up of the low frequency optical conductivity with lattice shown in Fig. |5j The 



data points in both curves are fit by the simple two-parameter Drude form (3.2) 



3.1 The Drude Peak 

At low frequency, both the real and imaginary parts of the conductivity can be fit by the two- 
parameter Drude form 

a{ U ) = T — (3.2) 

with both the scattering time r and the overall amplitude K constants, independent of u. This is 
shown in Fig. [6j It can be checked numerically that the overall amplitude K agrees (to about the 



1% level) with the coefficient of the pole (3.1) in the translationally invariant case. All interesting 
physics in this regime is therefore captured by the single parameter, r. We have varied the 
temperature and lattice spacing and found that this Drude form holds in all cases. Given the lack 
of well-defined quasi-particles in our holographic system, it seems surprising that the low-frequency 
behaviour of our system is governed so well by the exact Drude form. 

3.2 DC Resistivity 

The resolution of the u> = delta-function leaves behind a well-defined DC resistivity, p = (Kr)~ l . 
The Drude amplitude K is essentially independent of temperature T and all temperature de- 
pendence in the resistivity p(T) is inherited from r. The results depend strongly on the lattice 
wavenumber ko and are shown on the left hand side of Fig. [7j 

To make sense of this complicated plot, we review some recent work in the literature. Since 
the near horizon geometry of an extremal Reissner-Nordstrom AdS black hole is AdS^ x R 2 , the 
dual theory is said to be "locally critical" in the sense that it is invariant under rescalings of time, 
with no rescaling of space. Hartnoll and Hofman [12] have recently studied the DC conductivity 
in a locally critical theory. They showed that the DC conductivity can be extracted from the two 
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Figure 7: The left panel shows the DC resistivity plotted as a function of temperature for various 



lattice spacings. On the right hand side we factor out the scaling (3.3 3.4) and re-plot the same 



data on a log scale. The lines denote a fit to the data including polynomial corrections to the 
leading low temperature behavior. Both plots arise from a background with p = 1.4 and the 
lattice amplitude Aq = k^/2. The plots remain essentially unchanged for lattices of different 
amplitudes. 



point function of the charge density, evaluated at the lattice wavenumber. They then calculated 
this two point function by perturbing the Reissner-Nordstrom AdS black hole and found 

p oc T 2v - 1 (3.3) 

where^ 

v = \ ^ + 2{k/pf~^l + {k/pf (3.4) 

The exponent can be viewed as arising from the dimension A = v — \ of the operator dual to the 
charge density in the near horizon AdSi region, evaluated at the lattice wavenumber k. 

On the right hand side of Fig. [7] we plot pjT 2v ~ x for several values of the lattice wavenumber. 
As discussed earlier, if our scalar field has lattice wavenumber k , the charge density has lattice 



wavenumber 2ko, so we have set k = 2fc in (3.4). We have fit the data to po = T 2 "^ 1 ^ + a{T + 
a-iF 2 + g^T 3 ) and drawn the curves on the right hand side of Fig. [7j The fact that the curves all 
approach nonzero, but finite, constants at low temperature shows that our data confirms the low 



temperature scaling (3.3) with exponent (3.4) predicted in [12] . 

Note that as the temperature goes to zero, the dissipation goes to zero and the DC resistivity 
vanishes. Thus the DC conductivity becomes infinite, as expected for a perfect lattice with no 
dissipation. 



This is a manifestly scale invariant form of the exponent that was found in |12j and was first derived in a 
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Figure 8: The modulus |er| and argument argcr of the conductivity. The background for both 
plots has wavenumber ko = 2, amplitude Aq = 1.5, chemical potential \x = 1.4 and temperature 



T/pL — 0.115. The line on the right is a fit to the power law (3.5) 



3.3 Power-law Optical Conductivity 

The fit to the Drude peak works well for lo/T < 1. However, for lo/T > 1, the optical conductivity 
exhibits a power-law fall-off in a "mid-infrared" regime, before reverting to the continuum result. It 
is convenient to use lot as a dimensionless measure of the frequency in this region. This is because 
the power law behavior falls roughly in the range 2 < ur < 8 for all the lattices we have examined, 
even those with different temperature, lattice spacing and amplitude. (For the temperature we use 
in Fig. 5, tT = 2.22, so lot = 2 corresponds to lo/T = 0.9.) In Fig. [8] we have plotted \a\ and the 
phase angle over this range of frequencies. The data is very well fit by 

\a(fj)\ = ^ + C (3.5) 

Meanwhile the phase angle has only small dependence on lo, but varies between 65° and 80°, as ko 
varies from 1 to 3. The slight variation in the phase angle is enough so that the real and imaginary 
parts of the conductivity do not individually follow simple power laws over the range indicated in 
Fig. |} 

The exponent of the power law does not depend on the particular choices we have made for the 
parameters in our model. To illustrate this, and to make the power law more manifest, in Fig. [9] 
we show (|<r| — C) vs lot on a log-log plot. On the left we show three different choices for the lattice 
wavenumber fco- On the right, we show three different temperatures. The fact that the curves all 
form parallel straight lines for lot > 2 shows the power law fall-off with exponent —2/3 is robust. 
Since the offset C depends on k and T, in Fig. [9] we have subtracted a different constant for each 
curve. 



different context in 
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Figure 9: The magnitude of the optical conductivity with the offset removed, on a log-log plot. 
On the left, the plot has T//i = .115, and shows three different wavenumbers: diamonds denote 
k = 3, the squares denote k = 1, and the circles denote k = 2. On the right, the plot has 
ko = 2 and shows three different temperatures: the diamonds have T//i = .098, the circles have 
T/fj, = .115 and the squares have T/fi = .13. In both plots, Ao/ko = 3/4. The fact that the lines 



are parallel for ur > 2 shows that the fit to the power law (3.5) is robust 



3.4 Comparison to the Cuprates 

Our results for the optical conductivity bear a striking resemblance to the behaviour seen in the 
cuprates where a mid- infrared power-law contribution has long been observed [IHl HU [19] . At low 



frequencies, u/T < 1.5, the optical conductivity takes the Drude form (3.2) if one also allows for r 
to scale linearly with both T and u [19]. But for u/T > 1.5, a cross-over to power-law behaviour 
takes place, in which \a\ ~ with 7 f» 0.65 and the phase is roughly constant around 60° [TTJ . 



This is shown in Fig. 10 



There are differences between our behaviour and that of the cuprates. Most notably, the DC 
resistivity that we obtain is nothing like the robust linear behavior characteristic of the strange 
metal regime; instead we find power-law behaviour with an exponent that depends on the lattice 
wavenumber. Moreover, in the mid-infra red regime our optical conductivity requires a constant 



off-set C in (3.5). No such off-set is seen in [TT] . 

In some sense, this discrepancy makes the agreement of the phase angle somewhat more sur- 
prising. For strictly power-law conductivity (i.e. with C = as seen in the cuprates), causality 
and time-reversal invariance in the form a(u) = <y*(—u) relate the exponent of the power law to 
the phase [UJ. However, with C / 0, there is no reason for the phase and exponent to be related. 

There have been previous theoretical observations of power-law fall-off in the optical conduc- 
tivity. It is seen in models of charges moving in a periodic potential subject to dissipation [2D], but 
the exponent typically depends on the details of the model of dissipation. In the context of the 
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Figure 10: The optical conductivity of optimally doped Bi2Sr2Cao.92YomCu20s+s- This plot is 
taken from [111. 



cuprates, Anderson suggested a power-law fall-off cr(u>) ~ w -7 on the basis of a Luttinger liquid 
model, with 7 = 2/3 arising from a coupling to a gauge field [2*Tj. 

Perhaps more pertinent for the present discussion, the universal power law observed in the 
optical conductivity, together with aw/T scaling, was associated to an underlying quantum critical 
point in [llj . Of course the starting point of our holographic model is a strongly interacting critical 
point, albeit with a scale introduced by the finite density Moreover, for small temperatures, T /i, 
our model exhibits an emergent locally critical point, reflected by the near horizon AdS% x R 2 
regime. 

However, a second explanation was put forward in [22] where it was argued that the a ~ 
behavior with 7 rj 0.65 was a generic prediction of electrons interacting with a broad spectrum of 
bosons. Our holographic model is certainly not short of such bosonic modes and it is possible that 
these are responsible for our observed behavior. 

Finally, within the holographic framework, a <t{uj) ~ co~ 2 ^ 3 power-law was shown to arise from 
probe charged matter interacting with a strongly coupled soup with dynamical exponent z = 3 

ma- 

4 Future Directions 

The introduction of a gravitational lattice in the simplest holographic model of a conductor has 
allowed us to explore the low-frequency optical conductivity in these models. At very low frequen- 
cies, ct(u) follows a simple Drude form. However, for intermediate frequencies, |er(w)| has a power 
law fall off (with constant offset) and its phase is approximately constant. Remarkably, both the 
exponent of the power law and the phase are consistent with data taken on some cuprates and 
are robust against changing all parameters of our model. We do not have a deep understanding 
of why this is happening and it would clearly be of interest to find an analytic derivation of this 
result. 

To get more insight into this result, there are a few generalizations that should be investigated. 
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First, we have added the lattice in only one direction. It is natural to ask if the power or phase 
change if we do a more realistic calculation with a lattice in both spatial directions. We expect 
the answer will be no, but the calculation is considerably more difficult as it now requires solving 
nonlinear PDE's in three dimensions. Second, one can ask if the power and phase change if we 
consider a five dimensional bulk dual to a 3 + 1 dimensional system. This is currently under 
investigation. 

We have chosen to introduce a lattice by adding a periodic source for a neutral scalar operator. 
This induces a periodicity in the charge density, but the chemical potential \i is kept constant. 
One can alternatively model an "ionic" lattice" by making \i itself a periodic function of x. (This 
is the approach taken in [TJ1 12H ES]-) In this case, the dual gravitational description contains only 
Einstein-Maxwell theory. Work on this alternative lattice is underway and results will be reported 
elsewhere. Preliminary calculations indicate that the tu^ 2 ^ 3 scaling is seen for this lattice also. 

We have taken just the first steps toward introducing lattice effects into holographic discussions 
of condensed matter systems. There are many future directions. For example, it would be inter- 
esting to determine the phonon spectrum. Even though the lattice is fixed in the UV, the bulk 
bulk field it generates is dynamical and fluctuates. These fluctuations are governed by quasinormal 
modes. 

Another open problem is to add a charged scalar field to our gravitational model. Without the 
lattice, it is known that this theory has a second order phase transition at low temperature to a 
superconducting phase. It would be interesting to study the effect of the lattice on the transport 
properties of the superconductor. 

Finally, one can add a probe fermion to the bulk. In the translationally invariant case, this 
was used to find Fermi surfaces. One can now study more realistic Fermi surfaces in the presence 
of a lattice. 

A Appendix: Boundary Conditions 

In this Appendix we present the conditions on the perturbation at the horizon which implement 
ingoing wave boundary conditions. This can be easily achieved if one writes the background metric 
in terms of ingoing Eddington-Finkelstein coordinates, and require h a b to be a smooth tensor, at 
the horizon, in these coordinates. This induces the following near horizon behavior 

h tt (x,z) = (1 - z)"fa [U° t \x) + (1 - z)h { t l\x) + 0((1 - z)f] 

h tx {x, z) = (l- z)-^[hfj(x) + (1 - z)~hS(x) + 0{{\ - z)f] 

h tz (x, z) = (l- zY'-^hfXx) + (1 - z)h\l\x) + 0{{l - z)) 2 } 

h xx (x, z) = {\- z)-^m{x) + (1 - z)hQ{x) + 0((1 - z)f] 

h xz (x,z) = (1 - z)- l -^[hW{x) + (1 - z)h<8(x) + 0{{l - z)f] 

h„(x,z) = (1 - z)- 2 -^[h^(x) + (1 - z)h<g(x) + 0((1 - z)) 2 } . (A.l) 

h yy (x,z) = (1 - z)-™[h$(x) + (1 - z)hg\x) + - z)f] 

b t (x, z) = (1 - zy^lb^ix) + (1 - + 0((1 - z)) 2 } 

b x {x, z) = (1 - z)-^m[b^\x) + (1 - z)b^(x) + 0((1 - z)) 2 } 
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b z {x, z) = {\- z)- l -M-)\bf{x) + (1 - z)b^\x) + 0{{l - z)) 2 } 

fj(x,z) = (1 - Z y^[fl i0) {x) + (1 - Z)fj {1 \x) + 0((1 - 2)) 2 ] 



where -P(l) is given in (2.9). The boundary conditions are then fixed by inputing the above 
expansion into the equations of motion, and solving them in a (1 — z) expansion off the horizon. 
This expansion will fix some relations amongst the ^ coefficients. For instance, at the horizon 
regularity demands: 

= (A-2) 

The remaining regularity conditions express how the W coefficients relate with some of the ^ 
coefficients and their tangential derivatives. These are far too cumbersome to be presented here. 
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